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1. Background 

The phase structure of QCD at finite density and temperature is of fundamental importance, 
and can be studied both experimentally and theoretically. Nevertheless, progress has been slow, in 
part because of the sign problem, which afflicts both lattice simulations [ 1 ] and phenomenologi¬ 
cal models [2], The problem in QCD is due to the fact that the fermion determinant is complex 
for typical gauge field configurations when the quark chemical potential ji is nonzero. In [3], we 
have shown that the consideration of complex saddle points provides a conceptually cohesive phe¬ 
nomenological model of QCD at finite T and fj,. Moreover, we have identified a new property of 
QCD at finite density, the occurrence of a disorder line, that may have observable consequences 
in experiment and/or lattice simulation. Some feature associated with the disorder line differen¬ 
tiate strongly between different phenomenological models, and may thus have an impact on our 
understanding of confinement. 

We consider an SU(N) gauge theory coupled to fermions in the fundamental representation. It 
is well-known that the Euclidean Dirac operator has complex eigenvalues when a nonzero chem¬ 
ical potential is introduced. This can be understood as an explicit breaking of charge conjugation 
symmetry The log of the fermion determinant, log det (/r, A), which is a function of the quark 
chemical potential /./ and the gauge field A, can be formally expanded as a sum over Wilson loops 
with real coefficients. For a gauge theory at finite temperature, the sum includes Wilson loops that 
wind nontrivially around the Euclidean timelike direction; Polyakov loops are examples of such 
loops. At jU = 0, every Wilson loop Tr pW appearing in the expression for the fermion determinant 
is combined with its conjugate Ti pW T to give a real contribution to path integral weighting. More 
formally, charge conjugation acts on matrix-valued Hermitian gauge fields as 

^ '■ Afi —>• — A^, (1.1) 

where the overall minus sign is familiar from QED, and the transpose interchanges particle and 
antiparticle, e.g., VT and W in SU (2). This transformation law in turn implies that A' exchanges 
the Wilson loop and the conjugate loop as shown in Fig. 1(a), so unbroken charge symmetry im¬ 
plies a real fermion determinant. When /./ / 0, Wilson loops with nontrivial winding number n 
in the X 4 direction receive a weight e n ^ while the conjugate loop is weighted by e~ n ^ as illus¬ 
trated in Fig. 1(b), and thus invariance under is explicitly broken. However, there is a related 
antilinear symmetry which is unbroken: Tr pW transforms into itself under the combined action of 
ffJU, where JU is the fundamental antilinear operation of complex conjugation. Thus the theory 
is invariant under JU even in the case 0. Note that JU itself is also a symmetry of the theory 
when jU = 0 , but it is explicitly broken when /.( 7 ^ 0 for the same reason as r S ; . symmetry is an 
example of a generalized 2? (parity-time) symmetry transformation [4]; theories with such sym¬ 
metries form special class among theories with sign problems. Even though fo UfU is a trivial trans¬ 
formation for each Wilson loop, it implies the well-known relation det (—[X,A^) = det (fl.AU for 
Hermitian A^, because acting A' on the fermion determinant is equivalent to changing jU to —/r as 
one can see from Fig. 1, and acting JU on the determinant simply gives the complex conjugation. 
The advantage of using is that it is more general, leading to more insight into the sign problem 
and applying to bosons as well as to fermions. 
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Figure 1: The Wilson loop trVV and the Polyakov loop tr/ J under the charge conjugation A and the complex 
conjugation JP. 


For phenomenological models, the existence of fff symmetry leads naturally to the consid¬ 
eration of complex but c ff ,ff -symmetric saddle points. Typically, such models require the mini¬ 
mization of some effective action F or effective potential Veff as a function of some set of fields. 
We will consider models with effective potentials that are class functions of the Polyakov loop P, 
depending only on the set of eigenvalues of P. 'ff -ff symmetry will map any saddle-point con¬ 
figuration aJj 1 * into another saddle point given by = —/ij, 1; ' with a corresponding connection 
between the actions of the two configurations: S (2} = S' ]> *. However, some field configurations 
are themselves ^JF-symmetric in that —Ajj is equivalent to A Ll under a gauge transformation. If a 
saddle point is fff symmetric, then its action and effective potential are necessarily real. A quick 
direct proof can be given: For such a field configuration, it is easy to prove that every Wilson loop 
is real and thus det (ju.A/j) is real and positive for a r ff .fff -symmetric field configuration. If a single 
< A ? JF-symmetric saddle point dominates the effective potential, then the sign problem is solved, at 
least for a particular phenomenological model. Such ‘if JF-symmetric saddle points have been seen 
before in finite density calculations [5], 

Let us consider the Polyakov loop P, a special kind of Wilson loop, associated with some 
particular field configuration that is 'ff.ff -symmetric. We can transform to Polyakov gauge where 
A 4 is diagonal and time-independent, and work with the eigenvalues Oj defined by 


P (x) = diag 


■ ■ 


e ‘ Q n(x) 


(1.2) 


where the Of s are here complex but satisfy £ ; Oj = 0. Because we are primarily interested in 
constant saddle points, we suppress the spatial dependence hereafter. Invariance under '£-'Af means 
that the set |— 6 j | is equivalent to the { Oj } although the eigenvalues themselves may permute. In 
the case of SU (3), we may write this set uniquely as 


{6 - iy ,-6 -iy, 2 iy}. 


(1-3) 


This parametrizes the set of , iff -symmetric SU( 3) Polyakov loops. Recalling the form of the 
diagonal Gell-Mann matrices A 3 and As, we see that in Polyakov gauge Of 0 corresponds to a 
nonzero real value for A%, while y f 0 corresponds to a purely imaginary value for A\. Thus a 
%'fff -symmetric saddle point requires analytic continuation of A\ along the imaginary axis. Notice 
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that both 

TipP = e ¥ 2 cos 6 +e~ 2v 


and 


TrfP 1 = e~'< , 2 co sO + e 2 ^ 


(1.4) 


(1-5) 


are real, but they are equal only if y/ = 0. In the usual interpretation of the Polyakov loop expec¬ 
tation value, this implies that the free energy change associated with the insertion of a fermion is 
different from the free energy change associated with its antiparticle. It is easy to check that the 
trace of all powers of P or P' are all real, and thus all group characters are real as well. 

The existence of complex -symmetric saddle points provides a fundamental approach to 
non-Abelian gauge theories that is similar to the heuristic introduction of color chemical potentials, 
and naturally ensures the system has zero color charge, i.e. , all three charges contribute equally [6]. 
In the case of SU (3), extremization of the thermodynamic potential with respect to 6 leads to the 
requirement (n r ) — {«,,) = 0 where (/?,.) is red color density, including the contribution of gluons. 
Similarly, extremization of the thermodynamic potential with respect to i// leads in,) + {«.,,) — 
2 {rib) = 0- Taken together, these two relations imply that ( n r ) = ( n g ) = {rib). 

We demand that any saddle point solution be stable to constant, real changes in the Polyakov 
loop eigenvalues, corresponding for 51/(3) to constant real changes in A\ and A\ j. Consider the 
{N — 1) x (A — 1) matrix M a b, defined in Polyakov gauge as 


M ab = g 2 


d 2 V eG 
dA“dA b ’ 


( 1 - 6 ) 


At very high temperatures and densities, the eigenvalues of this mass matrix give the usual Debye 
screening masses. The stability criterion is that the eigenvalues of M must have positive real parts. 
At Tf./C-symmetric saddle points, the eigenvalues will be either real or part of a complex conjugate 
pair. In the case of SU (3), the matrix M may also be written in terms of derivatives with respect to 
6 and y/ as 



1 'TV’efl i (TCff \ 
4 d9 2 4vT dddy t 
i d 2 Vtff -1 'TV’eft / 
4f3 dOdy 12 dy 2 / 


(1.7) 


This stability criterion generalizes the stability criterion used previously for color chemical poten¬ 
tials, which was d 2 V e ff/dy r 2 < 0. Crucially, the mass matrix M uh is invariant under M* = O 3 M 03 , 
which is itself a generalized & (parity-time) symmetry transformation [4]. It is easy to see that 
this relation implies that M a b has either two real eigenvalues or a complex eigenvalue pair. In either 
case, the real part of the eigenvalues must be positive for stability. In the case where there are two 
real eigenvalues, we will denote by fq and K 2 the two positive numbers such that K 2 and fc| are the 
eigenvalues of the mass matrix M a b. If M a i, has two complex eigenvalues, we define two positive 
real numbers Kr and Kj such that ( K/ { ± iK/) 2 are the conjugate eigenvalues of M ( ,b. The border 
separating the region iq /0 from the region K 7 = 0 is known as the disorder line. In this case, it 
separates the region where the color density correlation function decays exponentially in the usual 
way from the region where a sinusoidal modulation is imposed on that decay. 
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2. Models 


We now consider a class of phenomenological models that combines the one-loop result with 
the effects of confinement for the case of SU(3) gauge bosons and two flavors of quarks at finite 
temperature and density. The model is described by an effective potential which is the sum of three 
terms: 

VeffOP) = V g (P)+Vf(P)+V d (P). (2.1) 


The potential term V g ( P ) is the one-loop effective potential for gluons. The potential term Vj (P) 
contains all quark effects, including the one-loop expression. The potential term V,i IP) represents 
confinement effects. We will consider three different forms for Vf ( P ) and two different forms for 
Vd (P) for a total of six different models. The formulas and parameters we use for these models can 
be found in [3]. 

The potential term V ( i(P) is taken to respect center symmetry and acts to favor the confined 
phase at low temperature and density [7, 8]. The gauge contribution V g (P) favors the deconfined 
phase, and in the pure gauge theory (Nf = 0) the deconfinement transition arises out of the compe¬ 
tition between V g (P) and Vd(P). The parameters of V t j(P) are set to reproduce the deconfinement 
temperature of the pure gauge theory, known from lattice simulations to occur at Td ~ 270MeV. 
The specific forms used here are Model A and Model B of [7], which can be written as 


Vd(P) 

v d B (P) 


N 


Ed 

j,k— l 






M\T 2 ((27T-30) 2 -27yr 2 ) 

6 n 2 



1 

4 


{cost? — cosh(3t//j} 2 sin 2 6 


( 2 . 2 ) 

(2.3) 


where we have used the parametrization in Eq. (1.3). A 0j k = \ Qj — 0 k | are the adjoint Polyakov loop 
eigenvalues and Z ?2 is the second Bernoulli polynomial. The expression for Model A gives a simple 
quartic polynomial in the Polyakov loop eigenvalues for V g (P) + V t f (P) and thus can be thought 
of as a form of Landau-Ginsburg potential for the Polyakov loop eigenvalues. The parameter M,\ 
controls the location of the deconfinement transition in the pure gauge theory, and is set to 596 MeV. 
On the other hand, the form for Model B is motivated by Haar measure, representing a determinant 
term that tries to keep a space-time volume of order j 8 R? confined. In order to reproduce the correct 
deconfinement temperature for the pure gauge theory, R must be set to R = 1.0028 fm. At low 
temperatures, the potential Vd(P) dominates the pure gauge theory effective potential. The variable 
i// is zero, and Vd (P) is minimized when 6 = 2n/3. For this value of 0, the eigenvalues of P are 
uniformly spaced around the unit circle, respecting center symmetry, and Tr pP = Try P' = 0. As 
the temperature increases, V g (P) becomes relevant, and gives rise to the deconfined phase where 
center symmetry is spontaneously broken. The addition of light fundamental quarks via V/(P) 
explicitly breaks center symmetry. For all nonzero temperatures, center symmetry is broken and 
(Tr pP) / 0. However, a remnant of the deconfinement transition remains in the form of a rapid 
crossover from smaller value of Tiy P to larger ones as T and /./ are varied. 

Although V c f and V c f appear to be very different, and are motivated in different ways, they are 
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actually closely related. The confining potentials V c i can also be written as 





V/ 



(2.4) 

(2.5) 


up to some irrelevant constant for the latter case. Using Tr^P = Tiy.-P''Try P : " — 1, it is easy to prove 
that minimizing either v£ or V - / yields a confining phase where Tr pP n = 0 for all n / Omod(V). 

We consider three different cases of quarks with mass m. The first is heavy quarks, with a 
fixed mass of 2 GeV. In this model, the quarks are essentially irrelevant for the deconfinement 
transition, which occurs at essentially the same temperature as if no quarks were present at all. 
The effect of spontaneous chiral symmetry breaking is not included, as it would only contribute a 
small amount to the quark mass. This case is in some sense the simplest, and perhaps would be 
the easiest for which to obtain reliable simulation results. The second case considered is massless 
quarks, where the fermion mass is set equal to zero by hand. This case cannot be easily simulated 
using lattice methods, because it ignores chiral symmetry breaking effects which do occur in lattice 
simulations. It is thus useful only for sufficiently large values of T and /./ such that chiral symmetry 
is essentially restored. Our most realistic treatment of quarks uses a Nambu-Jona Lasinio four- 
fermion interaction to model chiral symmetry breaking effects, so these models are of Polyakov- 
Nambu-Jona Lasinio (PNJL) type [9]. 


3. Disorder lines in phenomenological models of QCD 

We now show some results of the disorder lines, which are the borders separating the region 
Kj / 0 from the region kj = 0: see the last paragraph in Sec. 1. More results and discussions of the 
disorder lines and the Polyakov loops at the -symmetric saddle point can be found in [3]. 

Figure 2 shows contour lines for ly in the fi — T plane along with the region where Kj / 0 as 
well as the critical line for the PNJL model with Model A and Model B. In both models, the shaded 
regions where Ky / 0 cover a large portion of the phase diagram. While the critical line does not 
depend much on the models, there is a striking difference in the disorder line. The critical line lies 
completely within the region Ky / 0 for Model A while it appeases to be a smooth continuation of 
the critical line out of the critical end point for Model B. In both models, \j/ is very small and a jump 
in i// is visible as the critical line is crossed. Figure 3 shows similar plots as Fig. 2 but with contour 
lines for Kj in the fi — T plane. Kj jumps at the critical line only for Model A, because the critical 
line lies on the disorder line for Model B. Comparison of the two figures shows that the peak in ly 
occurs at a lower value of f than the peak in Kj, with the peak in i// occurring near (f = 200 MeV, 
T = 110 MeV) for Model A and {f = 250 MeV, T = 140 MeV) for Model B. 

The behavior of the disorder line for large T and f for Model A is known analytically [3]: 



This behavior is generic to Model A when T,fi 3> m, as one can see in Figs. 2-4. On the other 
hand, there is no disorder line in the limit r,/t> m for Model B. This is consistent with the fact 
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(a) Model A (b) Model B 

Figure 2: Contour plot of t ft in the p T plane for the PNJL model. The region where Kj / 0 is shaded. 
The critical line and its endpoint are also shown. 



(a) Model A (b) Model B 

Figure 3: Contour plot of Kj in the p T plane for the PNJL model. Contours are given in MeV, with a s 
set to one. 



(a) Massless quarks (b) Massive quarks (2 GeV) 

Figure 4: A comparison of the regions where Ki f 0. There is no disorder line for the case of massless 
quarks with Model B. High-T and low-T approximations are also shown for massless quarks. 
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that there is no disorder line in the entire phase diagram for massless quarks with Model B. This 
is the only case we have considered where there is no disorder line. Nevertheless, as shown in [3], 
y/ is nonzero even in this case, with a peak value near the same point as in the PNJL model with 
Model B. For heavy quarks, both models have the disorder line as can be seen in Fig. 4(b). Their 
shape is very similar for smaller values of /r, suggesting that some universal behavior occurs in this 
region. However, the behavior is very different in the region where both T and fi are becoming 
large. Model A shows a continuation of the disorder line that follows the behavior for massless 
quarks, while for Model B the disorder line covers a finite region in p — T space. The overall shape 
of the disorder line is similar to that found in the PNJL models, but of course shifted to a much 
larger value p. 

4. Conclusions 

As we have shown, the sign problem in QCD at finite density makes it very desirable to 
extend real fields into the complex plane. This extension is certainly necessary for steepest descents 
methods to yield correct results. Complex saddle points lead naturally to ( TrP) / (TrP^), a result 
that is much more difficult to obtain when fields are restricted to the real axis. The nature of these 
saddle points are restricted by '£-TP symmetry. The case of a single dominant saddle point is 
particularly tractable in theoretical analysis. In the class of models we have examined, the saddle 
point is not far from the real axis, as indicated by the small values of yf and corresponding small 
differences between (TrP) and (TrP'). This is good news for lattice simulation efforts, as it 
suggests only a modest excursion into the complex plane is needed. The small value of (// also 
indicates a small difference for thermodynamic quantities such as pressure and internal energy 
between our work and previous work on phenomenological models where only real fields were 
used. For all six cases studied here, the maximum value of y/ occurs in the region where quark 
degrees of freedom are “turning on,” as indicated by crossover or critical behavior. In our previous 
work on Model A for massless quarks [3], we were able to show analytically how y/ / 0 can arise 
from the interplay of confinement and deconfinement when Jt / 0. For the two PNJL models, it 
is striking that the largest values of y/ occur near the critical end point. These predictions can be 
checked in lattice simulations by the direct measurement of (TrP) and (TrP') once sufficiently 
effective simulation algorithms are developed. 

In all six cases studied, y/ / 0 leads to two different eigenvalues for the 4 4 mass matrix. In 
five of the six cases studied, a disorder line is found. This disorder line marks the boundary of the 
region where the real parts of the mass matrix eigenvalues become degenerate as the eigenvalues 
form a complex conjugate pair. In the PNJL models, the disorder line is closely associated with 
the critical line. Inside the region bounded by the disorder line, the complex conjugate pairs gives 
rise to color charge density oscillations. Patel has developed a scenario in which such oscillations 
might be observed experimentally [10]. Our results indicate that the oscillations may have too large 
a wavelength to be directly observable in experiment, although estimates based on phenomenolog¬ 
ical models should be applied cautiously. The mass matrix eigenvalues are in principle accessible 
in lattice simulations via the measurement of Polyakov loop correlation functions. A direct deter¬ 
mination of Kj may be difficult, but the disorder line itself could be determined from the merging 
of the values of Re( K\ ) with ReijCi ). 
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While the behavior of the Polyakov loop and the chiral condensate, as determined by lattice 
simulations, do not strongly differentiate between the two confining potential terms. Model A and 
Model B, the corresponding two-point correlation functions do. The most physically relevant case 
of PNJL models show both common features as well as clear differences in the behavior of the 
disorder line between Model A and Model B. In both cases, the maximum value of Kj occurs 
slightly above and to the left of the critical end point in the — T plane, in the vicinity of the 
region where the ratio TiyP' /TiyP is largest. In Model A, the critical line is contained within the 
boundary of the disorder line, but in Model B the disorder line appeals to come out of the critical 
end point as a continuation of the critical line, a common behavior for disorder lines. Furthermore, 
in Model A the disorder line continues diagonally in the ji — T plane for large /./ and T, but for 
Model B, the line bends over into the critical line. With Model A there is thus a possibility that 
the effects of the disorder line might be visible in the results of the Compressed Baryonic Matter 
(CBM) experiment at FAIR. The disorder line also strongly differentiates between Model A and 
Model B in the case of heavy quarks, so lattice simulations of either light or heavy quarks that can 
locate the disorder line have the potential to discriminate between the two models. 
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